function res = wishartln(V, W, v)

D = size(W,1);
denomln = wishartln_denom(W, v);
numerln = 0.5*(v-D-1)*log(det(V)) - 0.5*trace(inv(W)*V);

res = numerln + denomln;
